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Abstract. We propose a toy model with reversible mode coupling mechanism 
and with trivial Hamiltonian (and hence trivial statics). The model can be 
analyzed exactly without relying upon uncontrolled approximation such as the 
factorization approximation employed in the current MCT. We show that the 
model exhibits a kinetically driven transition from an ergodic phase to nonergodic 
phase. The nonergodic state is the nonequilibrium stationary solution of the 
Fokkcr-Planck equation for the distribution function of the model. 



1. Introduction 

First-principles understanding on the rich dynamic phenomena and the nature of the 
liquid-glass transition still remains as a challenging aim Q. As the only existing 
first-principle theory, the mode coupling theory (MCT) of supercooled liquids and 
the glass transition enjoyed considerable success in describing the dynamics of weakly 
supercooled regime of liquids Notwithstanding this surprising success, there are 
the following several unresolved issues concerning the basis of MCT: (a) A crucial 
ingredient of MCT is the factorization approximation which replaces the four-body 
time correlation functions by the product of two-body time correlation functions. This 
approximation is completely uncontrolled and its region of validity is a priori unknown, 
(b) The idealized MCT predicts a sharp dynamic transition to a nonergodic state at 
a certain temperature. But MCT does not provide any information on the nature of 
this nonergodic state, (c) The physical picture of the so called hopping processes in 
an extended version of MCT is still lacking. 

In recent years, possible deep connection between the structural glass and a class 
of spin glass models has been pointed out ||. In particular, the Langevin dynamics 
of the spherical p-spin model can be analyzed exactly in the thermodynamic limit 
due to the mean field nature of the model (i.e., full connectivity of the spins) Q |[. 
This analysis shows that the dynamic equation for the spin correlation function in 
equilibrium for p — 3 has the same form as in the Lcuthcusscr's schematic mode 
coupling equation for the density correlator || . The sharp dynamic transition observed 
in this class of models are driven by the dissipative nonlincarity in the equation of 
motion which originates from the nonlinear Hamiltonian J7j. In contrast to this, the 
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glassy behavior in the above-mentioned MCT (as well as in our model given below) 
is driven by the reversible nonlinearity [g which is dynamically generated and hence 
a non-trivial Hamiltonian is not necessary. 

With these situations, we thought that it is important to develop a toy model 
with following three ingredients: 

• reversible mode coupling mechanism 

• trivial statics 

• mean-field type so that the model can be exactly solvable. 

We have proposed such a toy model in a recent publication j|] . Here we further analyze 
the model. The model yields the self-consistent equations for the relevant correlation 
functions of the type familiar in the mode coupling theories of supercooled liquid and 
glass transition, where the strength of the hopping processes can be readily tuned. In 
the sense that the glassy behavior in this toy model is induced by the kinetics of the 
reversible mode coupling mechanism, our model is similar in spirit to the kinetically 
constrained models, the theme of the present workshop. 



2. Model 



Our model is defined as the following Langevin equations for the iV-component density 
variable ai(t) with i — 1,2, • • • , N and the M-component velocity variable b a with 
a = i, 2, • • • , M. Here and after we will use Roman indices for the component of a and 
Greek for that of b. 

d{ K{ a b(x -\- Jij a ajb a (1) 

K = -jb a - ui 2 K ja aj j==Jij a (uj 2 aiaj - TSy) + f a (2) 

v N 



< f a (t) > = 0, < f a (t)Mt') >= 2 1 TS a p5{t - 1') (3) 

where the summation is implied for repeated indices. Here 7 is the decay rate of the 
velocity b a and ui gives the j-independent frequency of oscillation of the density dj. 
The thermal noise f a (t) are independent gaussian random variables with zero mean 
and variance 27T, T being the temperature of the heat bath with which the system 
has a thermal contact. The choice of this variance guarantees the proper equilibration 
of the system. The N x M matrix Ki a plays an important role in the model and for 
later purpose we impose the (one-sided) orthogonality 

K ia Kif3 = S a /3, K ia Kj a 7^ 5ij (4) 

where the last equation is due to the inequality M < N. For M — N we can impose an 
additional condition Ki a = 6i a and hence trivially Ki a Kj a — Sij. We also note that 
Ki a governs linearized reversible dynamics of the model with the dynamical matrix 
ft given by fiy = uj 2 Ki a Kj a . The reversible nonlinear mode coupling terms are the 
ones involving the mode coupling coefficients Jij a which are chosen to be quenched 
(time-independent) gaussian random variables with the following properties: 

— J -0 



i] a 

2 f 

(SikSji + 5ii5j k )5af3 + Ki(3(K ka 5.ji + Ki a Sj k ) 
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JijaJkl/3 — „ 



+ Kjp(K ka 8 a + Ki a 5ik) 



(5) 
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where 7TtJ denotes average over the J's. Note that there is no thermal noise which acts 
directly on the density variable in (|l|) . This is because the model is constructed so as 
to mimic the dynamics of fluid. Equation (|l|) is analogous to the equation of continuity 
of fluid and Eq. @ is like the equation of motion where the right hand side is like the 
force acting on a fluid element. In constructing this model, we were motivated by the 
works 1 10, ll| in which random coupling models involving an infinite component order 
parameter have been shown to be exactly analyzed by mean-field-type concepts. We 
will thus eventually take N and M infinite with the ratio S* = M/N kept finite. 

One can derive from the Langevin equations (l)-(3) the corresponding Fokker- 
Planck equation for the probability distribution function D({a}, {&}, t) for our variable 
set denoted as {a}, {b} as follows 

d t D({a},{b},t) = LD({a},{b},t) (6) 

where the Fokker-Planck operator is given by L = Lq + L\ + Lmc 
with 

d / d , \ s / d , d 



L ° S db~J [ T db- a + K ) ' Ll = **« {-da- ba ' db a 
Lmc = —r^Jija { --^—ua.jb a + -^-w(o; 2 a l a J - T5«) 

VN V da i db a 



w 2 a- 



(7) 



It is then easy to show that the equilibrium stationary distribution (i.e., LD e (a, b) = 0) 
is given by 

D e ({a},{b}) = cst.e 2^,=i *r°i 2^ a=1 » J . ( 8 ) 
where est. is the normalization factor. 



3. Analysis and discussion 

For the subsequent analysis it is most convenient to introduce the following generating 
functional 

Z{h a ,h a ,h b ,h b }= [ d{a} ( d{b} [ d{a} [ d{b}e [ f dt{h >i +h i a i +h « ba+hbJa) e S (9) 



where the integrals are the functional integrals over the variable sets {a}, {a}, {&}, {b} 
and the h's and the /i's the conjugate source fields. The action S was decomposed 
into two parts Sq and Si which take the form 

S = Jdt^idi^di-Kiab^+iba^ba+jba+^K.aai-fa)^^) (10) 
Si = JjkaXjka (11) 

Xj ka = -^= J dt ^-idja k b Q + ib a w 2 aja k ^ (t) (12) 

where we have dropped the term TSij coming from eq. (||) since this term is negligible 
in the limit of infinite M and N. The functional determinant associated with the 
Langevin equations (0)-((^) which should appear in the integrand of the generating 
functional Z was equated to unity assuming the Ito calculus p2[ . The various 
correlation functions and response functions are obtained by taking various functional 
derivatives of In Z{h a , h a , h b , h b } with respect to h's and h's and setting them equal to 
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zero in the end in the standard way, where Z is the generating functional Z, averaged 
over the /'s and the J's. 

We now note that the replacements ia,j —> (uj 2 /T)aj, b a — > b a /T in Xj^ a leads 
to Xjk a = 0. Hence we can rewrite X^ a also as 



Xjka = Xjk a = -j= J dt ^-iajakb a + iw 2 & a <Zjafcj (t) 



(13) 



where ia, = i&i + (u> 2 /T)ai and ib a = ib a + b a /T. 

We now obtain for this toy model the equilibrium correlation functions defined 

as 

C a (t - t') = i < a 3 {t)a {t') >, C ab {t - t') = ^K ja < aj(t)b a {t') >, 



C ba {t - t') = —K ja < b a {t) aj (t') >, C b (t -t') = — < b a (t)b a (t r ) >, (14) 
Cf (t - = -^K la K ]a < a t (t)a 3 (t') > 



^K ja < b a {t)a 3 (t') >, C b (t ~t') = ^ 
1 

AT 

It turns out that we need to have the last correlation function to close the self- 
consistent set of equations for the correlators when M < N. Note that for the case 
M = N, if K ia — Si a is imposed, then (t — t') = C a (t — t'). The corresponding- 
response functions can be defined as 

G a (t - f) = 1 < o J -(t)ifl i (t / ) >, G afc (t ~ = -^ifja < a.j{t)ib a {t') >, 

G ba (t - f ) = — Jf jQ < 6 a (t)ia 3 -(t / ) >, G 6 (* - f ) = < b a {t)ib a (t') >, (15) 

Gf (t - = ^K ia K ja < ai (t)ia 3 (t') > 

Since we have a Gaussian stationary solution, we get the fluctuation-dissipation 
relationships (FDR) of the form 1 1 3 1 

G«(t - <') = - t')^C a (t - t'), G ab (t - = -0(t - O^C o6 (t - 0. 

where 0(t) is the unit step function: 0{t) = 1 for t > and otherwise. Note 
that this form of the FDR is rather unusual since the FDR usually takes the form 

G(t) = -e(t)d t c(t)/T. 

Another useful property arising from the causality and the above FDR is the 
following property 

< A{t)X{t') >=< X(t)A(t') >= for t > t' (17) 

for A(t) = (a{t),b(t)) and an arbitrary function X(t) = X(a(t),b(t),a(t),b(t)). 

We now take averages of Z over the thermal noise f a and the quenched random 
coupling Jija- In so doing we use the following properties which hold for the gaussian 
random variables: 

e -if dtt a (t)f a (t)\ = e - 7 T J dth a (t) 2 

. (18) 
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Defining the actions So and Sj as 

e So=^ e S \ e 5j =e^ J , (19) 

we obtain 

S a = J dt |iai(di - K ia b a ){t) + ib a (b a + -yb a + uj 2 K ia ai){t) - tTS^(*)| 

= J dt jifii ^ia, - TK la ib^j (t) + ib a (rib a + TKl-mi + jTib a ^ (t) J (20) 

where the last line is obtained using the property (|l7|). Now we have to deal with 

the interaction part Si = JjkaJimp Xjk a Xi m p jl. One can show that in the limit of 
M, N — > oo fluctuations can be neglected so that quantities like aj(t)a,j(t')/N etc. are 
replaced by C a (t,t'), etc. The interaction part Si then becomes gaussianized in the 
limit of M, N — » oo. The final expression for Si is then given by 

Si = I dt\iai(t)—^T, a a ® iOi(i) + K ia iai{t)TT, ab ® ib a (t) 



(22) 



+ K la ib a (t)^-T, ba <gi ioi(t) + ib a {t)TY> bb ® ib a (t) 1 (21) 
or J 

where X eg) a(i) = J_ <it' £(t — t')a(t') etc. Here the kernels E's are given by 
£ aa (t - t') = 5*^(C a (t - t')C b {t - t') + S*C ab (t - t')C ba (t - if)), 

2 4 

£ ab (i - t') = -25* 9 -^C a (t - t')C ba (t - t') 
H ba {t - t') = -25*^-C a (t - t')C ab (t - t'), 

Z bb (t-t<)^^^C a (t-tr 

Three kernels S oa , E ab , and £ ba comes from the nonlinear coupling term in the original 
Langevin equation (0) , and the kernel T, bb arises from the density nonlinearity in (g) . 
We also note that the correlator C^(t, if) is not involved in the E's. 

From the effective gaussian action S e f f = So + Si we can readily write down the 
following linearized Langevin equations for and b a 

di(t) = K ia b a (t) - S aQ ® Oi(i) - Jr to E«b ® &„(<) + /?(t) (23) 
6a(*) = -lK{t) - uj 2 K ia a,(t) - KiA ® ai (t) - E 6b ® 6 a (t) + /*(t) (24) 
where / a and / 6 are the effective thermal noises whose correlations are given by 

< ftWitf) > = ^[£~(tt / ) + E-(*'t)]<5<j 

1 



< fi(t)fa{t') > = K ia T[£ ab (tt') + -^Z ba (t't)} 



1 < 25 ) 
< /„(*)/*(* ) > = i^T[£ ba (tt') + — £ a6 (tt')] 
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Now we are ready to obtain a set of self-consistent equations for the five correlators 
from the linearized Langevin equations. By multiplying ( |23| ) by a,i(0)/N and (^) by 
K ia ai(0)/M and averaging over the effective thermal noise, we obtain 

C a {t) = 5*C ba (t)-V aa ®C a (t)-5*'E ab ®C ba (t) (26) 
Cba(t) = ~iC ba {t) - uj 2 C*{t) - E ha (8 Cf (t) - S fcb ® C ba (t) (27) 

where we used the causality requirements < ff (f)ai(O) >= and < /„(£)a»(0) >= 
0. Note that the correlator C% (t) appears in the equation for C ba (t). In order to 
obtain the equation for C^(t), we multiply ( |23| ) by KipKjpaj(0)/N and take thermal 
average. Then we obtain 

Ca (*) = C ba (t) - S Qa ® Cf (t) - S afc <g> C 6a (t) (28) 

Similarly by multiplying (p3| ) and (|24| ) by Kipbp(Qi)/M and b a (0)/M, respectively, and 
performing the thermal average we obtain the following equations for C ab (t) and C&(t) 

= C b (t) - S QQ ® C ab (t) - S Qb ® (29) 
C b (t) = - 1 C b {t)-w 2 C ab {t)-^ ba ®C ab {t)-^ bb ®C b {t) (30) 

The equations (p6|)-(p0|) constitute the self-consistent equations for the 5 correlators 
Ca(t), C ba (t), C£(t), C ab (t), and C b (t). This set of equations can be solved 
numerically with the initial conditions C a (0) = C*(0) = T/cj 2 , C ab (0) = C ba (0) = 0, 
and C b (0) = T. 

In analytic side, it is very convenient to work with the equations of the Laplace 
transformed correlation functions defined as C L (z) = J Q dte~ zt C(t). Performing the 
Laplace transformation of the self-consistent equations we obtain 



zC L a {z) 

zC« L (z) 

zC L ab (z) 
zC b L (z) 



— + 1 



-^ b (z))6*Ct(z)-^ a (z)Ct;(z) 
i- X*; b (z))Ct a (z)-(u 2 + Xt(z))C« L (z) 
^ b (z))Ct(z)-^ a (z)C« L (z) 



-3+1 



(l-^ b (z))Ct(z)-^ a (z)C^ b (z) 

T (c 2 + ^ a (z))CM ( 7 + ^M)C b L (z) 



From Q-pa), we obtain C£(z), C* L {z), and C^ a [z) in terms of E's as follows: 



T 



1 



u> 2 Z + ^ a {z) 



W a (l-E£ 6 (z)) 2 



(z + EJ (*))(* + 7 + £»(*)) + c 2 (l - (z))= 



T 



c 2 (l-S^(z)) 2 
z + 7 + £ bb (z) 
T{l-^ ab {z)) 



(31) 

(32) 

(33) 

(34) 
(35) 

(36) 
(37) 
(38) 

(39) 



[z + + 7 + £»(*)) + ^ 2 (1 - ^ b W) 2 

Here we have used the following symmetry relation 

^ a {z) = ~^ L ab {z) 

which follows from the definition of the kernels £ ab and £ ba , (l22|) , and C ab (t) = 
K ia < <H(t)b a {0) >= K ia < Oi(0)6 a (-t) >= < 6 o (*)ai(0) >= -C 6a (t). The 

first equality is due to the time translation invariance and the second one from the 
time reversal property of the velocity components. Note that for 6* = 1 the two 
correlators C^(z) and C^ L (z) become identical. 
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Similarly, from (p4[)-(p5D, we obtain 

a L (z\= ^X 1 ~ ^ab( Z )) IU1 

"'■ (z + Ei(z))( Z + 7 + S tt (z))+^(l-E^ 



C l (z) = T(z + ^ a (z)) 

b [ ) (z + EL (z))(z + 7 + + ^(1 - £^)) 2 

Now let us look at the behavior of the correlators for different values of S* . For 
S* = the only nonvanishing kernel is T, bb (z). Hence we obtain 

C L a {z) =~ ^ b (z)= 2 -^, (42) 



z 

,2 



z(z + 7) + (l + 2 5 2 T)u; 2 



(43) 



CbW = z(^ + 7) + (l + 2.g2T) W 2' ( 44 ) 

=<W = CTTFW (45) 

Here we point out that there appears to be a subtlety associated with the two 
limiting procesures: (A) first take 5* = M/N = before any calculation. (B) 
first calculate with S* > and then take the limit 5* — » 0+. The procesure (A) 
gives both C a (t) = C a {0) = T/uo 2 and Cf (t) = Cf (0) = T/w 2 . This is simply 
due to the fact that the {a} variables are time-independent since there is no velocity 
variable {b} that drives dynamics of {a}. However the results (|4^)-(^5]) were obtained 
by adopting the second limiting procedure (B). Here C a {t) is trivially nonergodic: 
C a (t) — C a (0) = T/uj 2 whereas and C^(t) exhibits a nontrivial nonergodic behavior: 
C*(t-> oo) = (T/Lu 2 )-2g 2 T/(l + 2g 2 T). The difference between these two procedures 
can be seen also by looking at ( p^ ) for C^ L (z). The terms except the first one on 
the right hand side is absent if the first limiting procedure (A) is adopted, whereas it 
remains finite in the second limiting procedure (B). 

For S* = 1 where M = N and K la = 5 la , C%(z) = C^ L (z) reproduces the 
equation derived in Jh|, apart from the wave number dependence. Note that if we 
put £aa( z ) = ^ab( z ) — by hand, (|3^ ) or (|37| ) gives a closed equation for C a (t) 
alone. This equation is nothing but the Leutheusser's schematic MC equation giving 
a dynamic transition from ergodic phase to nonergodic one. But in reality S QQ and 
E a (, can not be ignored and our numerical solution strongly indicates that the system 
remains ergodic for all temperatures due to the strong contribution of these so called 
hopping terms. Furthermore these hopping terms do not become self-consistently 
small as temperature is lowered. Therefore the density correlator does not show a 
continuous slowing down with lowering temperature. This result was striking to us 
since usually a mean-field-type theory, such as the dynamics of the spherical p-spin 
model in the limit of N — * oo, often gives a sharp dynamic transition. In fact, we were 
first constructing the toy model with M = N and we expected that the model designed 
to rigorously reproduce the idealized MCT exhibits such a dynamic transition. But to 
our surprise the dynamic transition was absent in the A^-component toy model. This 
aspect is a fundamental difference in the two kinds of mean-field-type theories with and 
without reversible mode coupling. The foremost example of the latter is the sherical 
p-spin model where the ergodic-to-nonergodic transition is driven by the dissipative 
nonlinearity which comes from the nonlinear random Hamiltonian. As demonstrated 




below, in order to have such a sharp transition in our toy model, we find it necessary 
to extend the original TV-component model to the model with M < N . Thus it 
is very difficult to understand the idealized MCT without relying upon uncontrolled 
approximation. It is also interesting to note that the ergodicity restoring process 
in our toy model (represented by the kernels S aa and £ h) has nothing to do with 
a thermally activated energy barrier crossing since the gaussian Hamiltonian in our 
model does not possess such a barrier. 

Our numerical solution for S* = 0.3 is shown in Figure 1 with various values of 
T. The other parameters were fixed as u> = 1, 7 = 1, and g = 1. As T is lowered, the 
relaxation exhibits a continuous slowing and it appears to be frozen at lowest T. One 
may ask whether this freezing reflects the presence of the genuine nonergodicity or it 
is merely apparent: the decaying will be observed if the observation time window is 
further extended. The question of the existence of nonergodicity is easily answered 
in the usual idealized MCT where one can easily solve the closed equation for the 
nonergodicity parameter to obtain the phase diagram. The situation is very different 
in our toy model. When we expand the correlators as C% (z) = f a / z + + z + ■ ■ ■ 
etc., we end up with a herarchically connected set of equations for all the f's, which 
can not be easily analyzed numerically. 

An analytic feature signifying the presence of the genuine nonergodic state can 
be seen by adiabatically eliminating the velocity components in the limit of large 
7 and obtaining the Fokker-Planck equation for the distribution function D({a},t) 
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containing only the {a} variables: 




(46) 




T 



(47) 



An important point is that the diffusion matrix Qij is singular for M < TV, i.e., 
det|Q| = @. The proof is simple. Define a N x N matrix M by My- = 
Mij =a for j < M, M,j = for j > M. Then we obtain in matrix notation 
Q = (T/7)M • M T (The superscript T denotes the transposed matrix). Then 
det|Q| = (T/7) Ar (det|M|) 2 = since det|M| = by construction. This implies 
that the Fokker-Planck equation ( [46|) can have nonequilibrium stationary solution 
other than the equilibrium one, D e ({a}) = est. exp(— Lu 2 a 2 /2T). This nonequilibrium 
stationary solutions are precisely the kind of nonergodic states found numerically in 
the present toy model. The general stationary solution ||l6| is given by 



where £j is the eigenvector of the diffusion matrix Qij with zero eigenvalue. If the 
function T{x) is a constant, then Dl{{o\) = D e ({a}) is the equilibrium distribution, 
otherwise it is a nonequilibrium stationary distribution. 

One instructive case for the nonequilibrium stationary solutions is that of 
g = 0. For this case, Qij becomes proportional to the dynamic matrix fi^: 
Qij — (T/j)Ki a Kj a — (T / ju> )Clij . By the same argument as above flij is singular 
as well. Note from @ that (z) = (T/lu 2 ) ■ (1 - S*)/z in the limit of z -> 0. 
The other correlators do not diverge at z = 0. Hence the model is nonergodic for 
< 8* < 1: the system is always driven into the nonergodic state in the linear case 
(g = 0). In this case the thermal noise alone is not enough to drive the system to 
the equilibrium state. This case is somewhat reminiscent of the ideal gas case or the 
collection of independent harmonic oscillators where the systems are trivially non- 
ergodic due to the absence of interactions. Only when the nonlinear reversible mode 
coupling is present, as T increases, the thermal noise can drive the system to the 
equilbrium state, hence making the system ergodic. The onset temperature at which 
the ergodicity is recovered is the dynamic transition temperature. 

In any event, further numerical and theoretical studies of possibile ergodic-to- 
nonergodic transitions for nontrivial case 5^0 are warranted. 

4. Summary 

We have constructed a dynamic mean-field-type model involving TV-component 
density and M-component velocity variables with reversible mode coupling and trivial 
Hamiltonian. The model is exactly solvable in the limit of N, M — * oo with keeping 
the ratio 6* = M/N finite. The model exhibits a sharp dynamic transition to a 
nonergodic state only in the range < 8* < 1. The nature of the nonergodic state 
can be understood in terms of the nonequilibrium stationary solution of the Fokker- 
Planck equation for the probability distribution for the density variable. It would be 
interesting to investigate the nonequilibrium aging behavior of the model. 



D L ({a})=F(Z j a j )e-^ 



(48) 
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